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ABSTRACT: In the minimal supersymmetric standard model, the conservation 
of i?-parity is phenomenologically desirable, but is ad hoc in the sense that it 
is not required for the internal consistency of the theory. However, if B — L is 
gauged at very high energies, -R-parity will be conserved automatically and exactly, 
provided only that all order parameters carry even integer values of 3(5 — L). 
We propose a minimal extension of the supersymmetric standard model in which 
.R-parity conservation arises naturally in this way. This approach predicts the 
existence of a very weakly coupled, neutral chiral supermultiplet of particles with 
electroweak-scale masses and lifetimes which may be cosmologically interesting. 
Neutrino masses arise via an intermediate-scale seesaw mechanism, and a solution 
to the problem is naturally incorporated. The apparent unification of gauge 
couplings at high energies is shown to be preserved in this approach. We also 
discuss a next-to-minimal extension, which predicts a pair of electroweak-scale 
chiral supermultiplets with electric charge 2. 



1. Introduction 



One of the successes of the Standard Model of particle physics is the automatic con- 
servation of baryon number (B) and total lepton number (L) at the renormalizable level. 
These conservation laws follow simply from the particle content and SU(3)(j x SU(2)l x 
U(l)y gauge invariance, and do not entail additional assumptions. 

The simplest supersymmetric extension of the Standard Model does not share this 
appealing feature, because the existence of scalar partners of the quarks and leptons allows 
for renormalizable violation of B and L. The most general renormalizable and SU(3)c x 
SU(2)i x ?7(l)y-invariant superpotential is given schematically by 

W = Wq + Wi + w 2 , 

W = [iH u H d + H u Qu + H d Qd + H d Le , 

Wi = udd , W 2 = \iLH u + QLd + LLe . 

[Here Q and L are chiral superfields for the 5 , L r (2)^-doublet quarks and leptons; u, d, 
e are chiral superfields for the SU(2) ^-singlet quarks and leptons, and H U: H d are the 
two SU (2) ^-doublet Higgs chiral superfields. It is possible to eliminate /j! by a suitable 
rotation among the superfields H d and L; but in general it is somewhat misleading to do 
so because in many extensions of minimal supersymmetry H d and L will not have the 
same quantum numbers, and the appropriate rotation cannot be performed.] The terms 
in Wo are just the supersymmetric versions of the usual standard model Yukawa couplings 
and Higgs mass, and they conserve B and L. However, W\ violates B by one unit and 
W 2 violates L by one unit. To prevent the proton from decaying within seconds or hours, 
either the couplings in W\ or those in W 2 (or both) must be extremely small. In this 
sense, the supersymmetric standard model appears to be less successful or at least less 
elegant than the Standard Model, since the observed conservation of B and L is no longer 
automatic, but requires some additional assumptions about the structure of the theory. 

The most common way to save the proton from the supersymmetric threat is to forbid 
all of the terms occurring in W\ and W 2 by imposing the discrete Z 2 symmetry [1,2] known 
as i?-parity or matter parity. The matter parity of each superfield may be defined as 

(matter parity) = (-lf( B - L ) . (1.1) 

Then multiplicative conservation of matter parity forbids all terms in W\ and W 2 , while 
allowing the phenomeno logically necessary ones in Wo. Equivalently, the .R-parity of any 
component field is defined by (— l) 3 ( B_Z/ )+ 2s j where s is the spin of the field. Since (— l) 2s 
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is of course conserved in any Lorentz-invariant interaction, matter parity conservation and 
-R-parity conservation are precisely equivalent. The description in terms of matter parity 
makes clear that there is nothing intrinsically ".R-symmetric" about this symmetry; in 
other words, it admits a formulation at the superfield level. Conversely, the description in 
terms of -R-parity is convenient in phenomenological discussions, because it happens that 
all Standard Model states have -R-parity +1, while all superpartners have -R-parity — 1. 
Conservation of -R-parity then immediately implies that superpartners can be produced 
only in pairs, and that the lightest supersymmetric particle (LSP) is absolutely stable. 

The minimal supersymmetric standard model (MSSM) with .R-parity conservation can 
provide a description of nature which is consistent with all known observations. However, 
the assumption of -R-parity conservation might appear to be ad hoc, since it is not required 
for the internal consistency of the theory. Alternative discrete symmetries have in fact been 
proposed (see for example [3,4]). Perhaps the simplest of these is the Z% discrete "baryon 
parity" of Ibahez and Ross [4] , which turns out to imply the falsifiable predictions that 
the proton is absolutely stable and there can be no neutron-antineutron oscillations even if 
there are isosinglet quark superfields near the TeV scale [5] . One might also entertain the 
possibility of small -R-parity violation, with intriguing phenomenological consequences (see 
for example [6-9]). However, if .R-parity is not exact, the LSP is unstable and so cannot be 
a candidate for the cold dark matter, unless its lifetime is of order the age of the universe. 

Fortunately, there is a particularly compelling scenario which does automatically pro- 
vide for exact -R-parity conservation due to a deeper principle. This is suggested imme- 
diately by (1.1), which shows that matter parity is simply a Z2 subgroup of B — L. If 
U(1)b~l 1S gauged at high energies, it will forbid each of the terms in W\ and W2 [10- 
13]. Of course, there is no massless gauge boson found in nature which couples to B — L, 
so U(1)b~l m ust be spontaneously broken. The question then becomes how to break 
B — L without also breaking matter parity. To guarantee that matter parity should re- 
main unbroken even after a gauged U(1)^-l * s broken, it is necessary and sufficient to 
require that all scalar vacuum expectation values (VEVs) or other order parameters carry 
3(B — L) charges which are even integers. Following the general arguments of Krauss 
and Wilczek [14], the gauged U(1)b~l symmetry breaks down to a Z2 subgroup which, 
in view of (1.1), is nothing other than matter parity. Unlike a global discrete symmetry, 
such a gauged discrete symmetry must be respected by Planck scale effects, and satisfies 
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discrete anomaly cancellation conditions [12,15,16]. Note that it is a contradiction in terms 
to speak of explicit .R-parity breaking in a supersymmetric model with gauged U(l)p_p; 
.R-parity will either be exactly conserved (if all order parameters carry only even integer 
values of 3(B — L)) or spontaneously broken [17] (if some order parameter carries an odd 
integer 3(B — L)). 

Of course, this scenario for the origin of matter parity is hardly mandatory, since 
it is technically natural to forbid the terms in W\ and Wi "by hand" as an unexplained 
assumption. However, it is worthwhile to take seriously the idea that -R-parity conservation 
is explainable, since we then obtain some quite non-trivial information about physics at 
very high energy scales. Not only do we gain an indication that the unbroken gauge group 
should contain U(X)b—Li but we also obtain information about how it should (and should 
not!) be broken. 

In ref. [13], the criteria for maintaining natural i?-parity conservation in models with 
gauged U(1)b~l were considered for various extended gauge groups. Consider, for exam- 
ple, the possibility of a natural explanation for i?-parity conservation in a supersymmetric 
SO (10) grand unified theory (GUT). Now, SO (10) contains U(1)b~l as a subgroup, so 
that i?-parity conservation is automatic before spontaneous symmetry breaking. How- 
ever, the smallest "safe" representations for a scalar which can break U(l)p_p without 
breaking the gauged matter parity subgroup in the process are 120, 126, 210 . . . and their 
conjugates. This is unfortunate, since experience has shown that it is quite difficult to 
build a successful supersymmetric GUT with such large representations. The alternative 
is to break 50(10) with an order parameter in a 16 representation, and that is what is 
usually done. However, the neutral component of a 16 carries 3(B — L) =3, so that 
the original automatic i?-parity conservation is forfeit. Indeed, renormalizable matter par- 
ity violation appears in the low energy superpotential from non-renormalizable operators 
of the form (1/M)(16) x 16 x 16 x 16. (Here and in the following M is some physi- 
cal cutoff scale, perhaps M = Mj>i anc k/y/87r.) As another example, one might consider 
an extension of the MSSM with a Pati-Salam gauge group SU(4)pg x SU(2)p x U(l)p. 
Again, i?-parity conservation is automatic before spontaneous symmetry breaking since 
SU(4)pg D SU(3)c x U(1)b—l- To avoid breaking matter parity in the process of break- 
ing SU(4)pg, it is necessary and sufficient that all order parameters have even SU(4)pg 
quadrality, since SU(4)pg quadrality = 3(B — L) [mod 4]. The smallest such "safe" rep- 
resentation for an order parameter which breaks U(1)b_l is the 10 of SU(4)pg. 

In ref. [18], the dynamical issues associated with automatic i?-parity conservation have 
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been considered in the case of left-right symmetric models. It was found that in a wide class 
of such models, .R-parity must be spontaneously broken because of the form of the scalar 
potential, although this can be evaded if non-renormalizable interactions are included. 

In this paper, we will consider instead a minimal extension of the MSSM in which the 
gauge group is extended by only U(1)b-l- Anomaly cancellation for U(1)b-l implies the 
existence of three neutrino chiral superfields v which carry B — L = 1 and are singlets of the 
standard model gauge group. A VEV for the scalar component of v would spontaneously 
break matter parity, so we will require it to be absent. While such a weak-scale VEV for v 
is not yet ruled out phenomenologically, we adopt for this paper the point of view that this 
is unacceptable, since we want to explore here only possibilities with exact and automatic 
-R-parity conservation. 

To obtain a realistic theory of neutrino masses, we may invoke the seesaw mechanism 
[19] by means of the superpotential 

W D y v H u Lv + — 5W . (1.2) 

Here S is a chiral superfield which must carry B — L = —2. Assuming that (S) is much 
larger than the electroweak scale m^, one finds that the lighter neutrino mass eigenstates 
have tiny masses ~ (y l/ {H u )) 2 /(ys{S)). Now, the role of S within this framework might 
be played by a composite field S = VV/M. However, this again cannot be consistent with 
our criteria for automatic -R-parity conservation, since then a VEV (S) ^ implies a VEV 
(17) ^ 0. Therefore, we prefer the possibility that S is a fundamental chiral superfield, 
so that the VEV (S) cannot break the matter parity subgroup of U(1)b-l- The field 
S must be accompanied by a field S in the conjugate representation, in order to cancel 
the anomalies and to allow spontaneous symmetry breaking in a nearly .D-flat direction. 
(Otherwise there would be catastrophically large supersymmetry-breaking D-terms, which 
would destabilize the electroweak scale.) For the models in this paper, the scale (S) is 
an intermediate one, roughly the geometric mean between the electroweak scale and the 
Planck scale. Assuming that the Yukawa couplings y v are of the same order as those of 
the charged leptons, one then expects light neutrino masses in the range relevant [20] to 
solar or atmospheric neutrino oscillations and hot dark matter. 

In section 2 of this paper we will propose a minimal extension of the supersymmetric 
standard model which successfully implements automatic and unbroken i?-parity conser- 
vation from gauged B — L, and discuss some of its implications. Section 3 contains some 
discussion of the subtleties associated with U(l) mixing in this model, and the effect of 
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intermediate scale thresholds on the sparticle spectrum. In section 4 we will discuss a next- 
to-minimal model [with gauge group SU(3)c x SU(2) L x SU{2) R x U(1)b-l\- Section 5 
contains some concluding remarks. 



2. A minimal model of automatic .R-parity conservation 

We consider a supersymmetric model with gauge group SU(3)c x SU(2) L x U(l)y x 
U{1)b-l- The MSSM chiral superfields (plus v) transform under this gauge group as three 
copies of 

Q ~ (3, 2, 1,1) d ~ (3, 1, J, -J) u ~ (3, 1, -§, - J) 

L ~ (1, 2, - J, -1) e ~ (1, 1, 1,1) v ~ (1, 1, 0, 1) 

and two Higgs doublets 

tf u ~(l,2,J,0) (1,2,-1,0). 

In order to break U{X)b-Li we introduce two chiral superfields 

S~ (1,1,0,-2) 5~(1,1,0,2). 

Because this gauge group contains two abelian factors, in general one must consider the 
possibility of arbitrary mixing between U{l)y and U(1)b-L- (Indeed, it would be quite 
unexpected if the gauge interactions were diagonal in the Y, B — L basis.) We will be able 
to postpone a discussion of this until the next section, however. 

Besides the interaction between S and v given in (1.2), S and S participate in a non- 
renormalizable superpotential interaction^ 

W = J^S 2 S 2 (2.1) 
and soft supersymmetry-breaking terms 

V so{t = m 2 \S\ 2 +m 2 \S\ 2 - (^ 2 S 2 + c.c.) . (2.2) 

(We use the same symbol for each chiral superfield and its scalar component.) The pa- 
rameters m 2 , m 2 and A should each be of order the electroweak scale in order not 



t We do not allow tree-level superpotential mass terms SS and H u Hd or their soft super- 
symmetry breaking counterparts, in sympathy with general results in superstring models. 



to upset the hierarchy. Note that the A/M term, while dimensionless, should nevertheless 
be treated as "soft" because of its tiny magnitude. Such terms should naturally arise in 
supergravity models, and this one will play a crucial role on several accounts, as we shall 
soon see. By a suitable phase rotation, we take A to be real and positive, while the phase 
of A can be arbitrary. The full scalar potential for the S and S degrees of freedom is given 
by the sum of Vg ft and 

%JSY = ^l^! 2 (|S| 2 + \S\ 2 ) + 9 -f (\S\ 2 - \S\ 2 f (2.3) 

where the latter term is the D-term contribution. The parameter gx is related to the 
gauge couplings of the U(l)s, and will be explicitly identified in the next section when we 
discuss the effects of U(l) mixing. 

A familiar method of inducing spontaneous gauge symmetry breaking in supersymmet- 
ric models is to arrange for the running soft (mass) 2 of the appropriate scalar to become 
negative at some scale. In the usual MSSM, this radiative symmetry breaking is achieved 
by means of a large top-quark Yukawa coupling which drives the Higgs (mass) 2 negative. 
In the model discussed here, the parameter m 2 obtains a negative radiative correction due 
to the Yukawa coupling yg in (1.2), and this has been exploited to obtain radiative sym- 
metry breaking in similar models [21-26]. However, in the present case one also obtains a 
large positive radiative correction to both m 2 and m 2 from £7(l)-gaugino loops. Indeed, 
the large (±2) B — L charges of S and S make it seem somewhat problematic to achieve 
radiative symmetry breaking in the traditional way; an examination of the renormalization 
group (RG) equations shows that very large (or numerous) Yukawa couplings yg seem to 
be required. Fortunately, it is not really necessary for m or m to be driven negative in 
this model, since the A term in (2.2) always favors spontaneous symmetry breaking. A 
non-trivial local minimum of the scalar potential will be obtained provided that 

A 2 -6|A| 2 (m 2 + m 2 ) >0. (2.4) 

This minimum will be global if 

A 2 - 8|A| 2 (m 2 + m 2 ) > . (2.5) 

These conditions can be satisfied either by driving m 2 negative, or simply by taking the 
free parameter A to be sufficiently large (while still roughly of order the electroweak scale), 
or perhaps by a combination of these effects. In any case, the minimum of the potential 
occurs along a nearly D-flat direction: 

(S) 2 w (S) 2 w ^2 ( A + \j A2 "6|A| 2 (m 2 + m 2 ^j (2.6) 
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with the deviation from .D-flatness given by 

(S) 2 -(S) 2 ^(m 2 -m 2 )/(2g 2 x ). (2.7) 

We see from (2.6) that the characteristic scale of B — L breaking is roughly a geometric 
mean between the electroweak scale and the Planck scale: (m^M) 1 / 2 ~ 10 10 GeV. Since 
only S and S obtain VEVs, it is clear that -R-parity conservation is automatic (and in 
fact unavoidable) in this model. Note that the minimum of the scalar potential is stable 
against v obtaining a VEV (and more generally against arbitrary perturbations of S, S and 
u); this can be understood from the fact that the scalar potential contains large positive 
semi-definite contributions |-Fj/| 2 = \yg(S)u\ 2 . 

After spontaneous symmetry breaking, a gauge boson and gaugino obtain masses 

M I = 2g x {S) (2.8) 

by eating the would-be Nambu-Goldstone boson degree of freedom Im[S — S]/V2. (With- 
out loss of generality, we take (S), (S) to be real and positive.) Here gx is the same 
coupling appearing in (2.3). In addition, one real scalar degree of freedom (given approxi- 
mately by Re[(5' — (S)) — (S — {S))]/y/2) and one Weyl fermion from S, S get masses Mj, 
forming a complete massive vector supermultiplet. There remains one light neutral chiral 
supermultiplet, given approximately by 

$^ [S- (S) + S- (S)]/V2, (2.9) 

whose components all obtain electroweak- scale masses. These degrees of freedom consist 
of a Weyl fermion ip (with .R-parity —1) of mass 

/o\2 

m ^ = 3|A| M (2 - 10) 



and two real scalar degrees of freedom a and b (of -R-parity +1) with squared masses 

m l = 4A ^ m 2 = 2A^-2(m 2 + m 2 ). (2.11) 



Let us pause to remark on several interesting features of this spectrum of electroweak- 
scale, neutral particles. First, note that in the limit A — > 0, m 2 vanishes and the scalar 
a becomes a Nambu-Goldstone boson. This corresponds to the spontaneous breaking of 
a continuous i?-symmetry of the superpotential, which is explicitly broken only by the A 
term. Fortunately, there is no reason for the parameter A to be small compared to the 
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electroweak scale; on the contrary, it is likely that A should be large in order to achieve 
the necessary condition (2.4) for spontaneous symmetry breaking, as we have already 
discussed. Perhaps a more plausible limit physically is A 3> |A| 2 (m 2 + m 2 ) (but still very 
roughly of order the electroweak scale), which leads to 

, ^ o AM A [A A [2 A , 

In general, the masses of the component fields of the supermultiplet $ satisfy the sum rule 

ml + ml - 2m 2 = m 2 + m 2 . (2.13) 

It is not difficult to show that the lightest member of the supermultiplet $ is always one 
of the scalars (a or b). 

An important byproduct of this symmetry breaking scenario follows from the existence 
of an allowed term in the non-renormalizable superpotential which is of the same order as 
(2.1): 

WD^H u H d SS. (2.14) 
After symmetry breaking, one obtains the usual fiHuH^ term of the MSSM, with 

H = \'(S) 2 /M (2.15) 

which is naturally of order the electroweak scale. The corresponding soft MSSM Higgs mass 
term (often denoted Bfi) is generated in the same way from the soft term corresponding to 
(2.14). This is a solution to the problem of generating an electroweak-scale \x term along 
the lines of [27]. 

The interaction (2.14) also plays another crucial role in this model; it allows a, 6, and 
ip to decay in a cosmologically timely fashion. These fields clearly have only tiny couplings 
to the particles of the MSSM. The most important such interactions actually follow from 
(2.14); one finds the coupling of $ D (a, 6, ip) to MSSM states 

WD WM* Hu H d . (2.16) 

The dimensionless coupling y^A' {S)/M is very roug hly of order (m^/M) 1 ^ if A' is of 
order unity, leading to decay widths of order V ~ m 2 ^/(87rM) ~ 1CP 15 GeV for the 
component fields of $ into MSSM states. Depending on kinematic and mixing angle 
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factors, one has^ two-body decays a, b — > Higgs + Higgs or a, 6 — > Higgsino + Higgsino 
and V' —> Higgs + Higgsino decays with lifetimes of order 10~ 9 seconds, give or take 
several orders of magnitude. (This should be compared to the notorious problem of the 
lifetime of an electroweak- scale gravitino, which can be estimated to be of order M/m\y 
times as long.) Thus it is unlikely that late decays of these particles could jeopardize the 
successful predictions of big-bang nucleosynthesis, although they might certainly have other 
interesting cosmological effects which should be carefully investigated. Note in particular 
that each decay of ip results in the production of one stable LSP. (It does not appear 
viable to allow ip itself to be the LSP, because its annihilation cross-section is so tiny that 
it would cause the universe to become matter dominated too early.) 

3. U(l) mixing and the sparticle spectrum 

The model described in the previous section contains two U(l) factors which can mix in 
an a priori arbitrary way. We chose to specify the charges of the chiral superfields in the Y, 
B — L basis, but this does not completely specify the gauge interactions of these fields. In 
fact, it would be rather surprising if the gauge interactions at high energies were diagonal 
in this basis. We will choose instead to use the "5'O(10)-inspired" basis given by U(1)r, 
U(1)b-Li with the B — L charges as before, and R = Y — (B — L)/2. [There would be no 
mixing in this basis if the gauge group were imbedded in e.g. unbroken SO (10).] One can 
always perform a rotation on the U(l) vector supermultiplets so that the kinetic terms are 
diagonal and canonically normalized for the U(1)r, U(1)^_i i gauge bosons and gauginos. 
Then the interactions with matter fields fa are specified by the covariant derivative 

D^fa = (fy + igfA% + igf- L A*- L )(t>i , 

gf = g R R l + g B - L , R y/3/8 (B - L), , 

9t L = 9b-lVV8 (B - L)i + g R ,B- L Ri ■ 

The charges R^ and (B — L)i are constants and are not renormalized. However, in general 
the couplings gs-L, g R , gB-L,R, and g Rt B-L all require counterterms and are renormalized 
[28]. The mixing couplings gB—L, R and g R ,,B~L cannot avoid counterterms unless the matter 
content is special, e.g. in complete multiplets of a non-abelian group containing at least one 
of the U(l) : s. It is therefore not consistent in general, and in particular in the model of the 
previous section, to set gs-L.R and g R .B—L equal to 0. At any particular renormalization 

t The components of $ can also decay into light (s) neutrino pairs, but these decays turn out 
not to be competitive because they are suppressed by the seesaw mixing angle squared. 
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scale one can perform a rotation on the vector superfield basis to set either g B - L , R or 
gn,B~L equal to [28]. This condition is not renormalization scale-invariant, however, so 
it is sometimes convenient to keep all four couplings as free parameters. 

In terms of these parameters, the coupling gx appearing in the previous section is 

g\ = (g R - y/3/2 g B . L<R ) 2 + {^/zj2 g B _ L - g R , B - L ) 2 . (3.1) 

At the scale of symmetry breaking, the surviving U{l)y gauge coupling is given by (in a 
GUT-like normalization) 



(3.3) 



g Y = a/572 (gagB-L - gB-L,Rg R ,B- L )/gx ■ (3-2) 

The one- loop RG equations for the gauge couplings are [t = \h(Q/Qq)]: 

gs-L 9b-l,r\ _ 1 ( gB-L 9b—l, R j ( b B _ L b R . B _ L 
dt \9r,b-l gR J 167T 2 \9r,b-l g R J \b R , B - L b R 

with 

b B - L = 9(gl-L + 9 2 r ,b-l) - zVGg B - L g R ,B-L , b R = 9{g 2 R + g B _ LjR ) - 2V§g R g B - L , R , 

b R , B - L = 9(g R g R ,B-L + 9b-l9b-l, r ) - V& (g B - L g R + 9b—l,r9r,b—l) ■ 
Using these equations, one finds that gy defined by (3.2) satisfies the one-loop RG equation 

d 1 33 o 

Jt' jY = m?r j y (34) 

(just as in the MSSM) both above and below M/, so that the condition for unification of 
gy with the SU(3)c and SU(2)l gauge couplings #3 and g<i is unaffected by mixing, up 
to two- loop and threshold effects. 

It therefore is sensible to impose a gauge coupling unification condition on all the 
couplings, as could follow from a superstring or a GUT model. At the unification scale tjj, 
one might therefore take 

93 = g2 = gR = gB-L = gu, 9b-l, R = 9 R ,b—l = . (3.5) 

We will assume these boundary conditions for the remainder of this section, although it 
cannot be overemphasized that alternative boundary conditions are certainly possible. At 
lower scales, one can then solve the one-loop RG equations analytically (for example by 
rotating to the multiplicatively renormalized basis), with the result 

9r = gB-L = gu («+ + K -)/ 2 » ( 3 - 6 ) 

g R , B -L = g B -L, R = gu («- - «+)/2 , (3.7) 



K± = 



1 + ^(9±V6) (tu-t) 



1-1/2 



(3.8) 



11 



(The first equality in each of (3.6) and (3.7) is due to a coincidental symmetry of the RG 
equations in this model.) On this "unification trajectory", the mixed couplings Qb-l,r and 
9r,b-l remain fairly small (< .04). 

Since the apparent unification [29] of gauge couplings observed at LEP can be main- 
tained in this model, it is sensible to explore features of the low-energy theory which follow 
from unified supergravity-inspired [30] boundary conditions. These boundary conditions 
include the supposition that at some scale Mjj > 2 x lO 1 ^ GeV the scalars in the theory 
have a common soft supersymmetry breaking (mass) 2 (denoted m 2 ,) and there is a common 
mass to -j/2 for each gaugino. One can then integrate the RG equations from Mjj down 
to the electroweak scale, and study the resulting low-energy theory. Here we will restrict 
ourselves to some brief comments regarding the impact of the extension of the MSSM of 
section 2, using the MSSM (with no new fields below Mjj) as a template. 

It is possible to show that the well-known gaugino mass unification prediction 

(M 3 /gl) = (M 2 /g 2 2 ) = (Mi/<&) = (m 1/2 /<#) (3.9) 

at low energies is precisely maintained by the one-loop RG equations of this model along 
the unification trajectory, provided that the gaugino masses are unmixed at the scale tjj 
in the R, B — L basis. [There is mixing induced among the gaugino mass parameters in 
the R, B — L basis by RG running, yet the surviving {7(1) y gaugino mass parameter does 
satisfy (3.9).] Therefore the predictions for chargino, neutralino, and gluino masses are 
essentially unaffected in the model of section 2, compared to the MSSM as a template. 
The condition (3.9) is modified by small two-loop corrections [31,32], of course. 

The predictions for masses of squarks and sleptons are affected in two ways. First, 
one has D-term contributions to scalar masses [22-26] due to the spontaneous breaking 
of the f7(l) symmetry. In the model of section 2, one finds that each scalar fa obtains a 
contribution to its (mass) 2 of 

Am? = (m 2 - m 2 ) [— X; + — Y^ 2 - g\_ L + g\ B _ L - g\_ L R (3.10) 



[9r9b-lm + g B -LgR,B-L])/gx\ 



where Xj = |YJ — |(J3 — L)j and m 2 — m 2 is evaluated at the scale Mj. For the unification 
trajectory, the term proportional to is non-vanishing but small, however it is important 
to keep in mind that it need not be so with more general boundary conditions on the gauge 
couplings. The corrections (3.10) should be added to the scalar masses at the intermediate 
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scale Mj and must be renormalized down to the electroweak scale; this turns out [26] to 
not affect the contributions proportional to X^, while inducing a quite small change in the 
term proportional to Y^. 

The presence of the additional U(l) gauge interactions above Mj also makes a con- 
tribution to scalar masses, because of terms in the RG equations due to gaugino loops. 
Evaluating these contributions for the unification trajectory of the RG equations (and 
taking into account all mixing effects) one finds the following approximate results for the 
slepton masses at the electroweak scale: 

m§ = [ml + .15m 2 , 2 - sin 2 d w m\ cos 2(3] - ^(m 2 - m 2 ) + .015m 2 /2 (3.11) 
m 2 ~ L = [ m 2 + .52m 2 / 2 + (sin 2 9 W - ^)m|cos2/3] + t^(™ 2 - m 2 ) + .045m 2 ^ 2 (3.12) 
TO ? = [to§ + .52m 2 ^ 2 + \m\ cos 2/3] + t^(™ 2 - m 2 ) + .045m 2 / / 2 (3.12) 

In each case the first set of terms in square brackets is the result for the template (MSSM) 
model. Next is the D-term associated with breaking of the £7(1) gauge group (neglecting 
the small contribution proportional to Y, and with to 2 — m 2 evaluated at Mj), and the 
last term is the additional contribution from U(l) gaugino loops above Mj. We have used 
here representative values Mjj = 2 x 10 16 GeV and Mj = 10 10 GeV. Similar equations 
can be written down for the squarks, although the relative effect is much larger for the 
sleptons. As long as we are assuming scalar mass unification, there is good motivation for 
the expectation that m > m , since m receives a negative RG contribution proportional 
to \ys\ 2 [cf. eq.(1.2)]. Therefore, the change (compared to the MSSM) in the difference 
between charged slepton masses, 

A(mj L - m\ R ) w (m 2 - m 2 )/5 + .03mJ /2 (3.13) 

is expected to be positive. This reinforces the expectation in the MSSM that should 
be heavier than e R ; depending on the relative magnitudes of m 2 , m 2 /2 , and m 2 - m 2 , the 
difference (3.13) could even be dramatic. 

4. An extension of the minimal model 

The model described in the section 2 is the simplest extension of the MSSM with 
gauged B — L breaking to matter parity at an intermediate scale. Its other successful 
features include a natural solution to the fx problem and a potentially successful theory 
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of neutrino masses. It is interesting to consider extensions of the minimal model with 
an enlarged gauge symmetry at the symmetry breaking scale. Note, however, that larger 
gauge groups seem to be somewhat disfavored for the following reasons. In order to have a 
viable seesaw mechanism, the order parameter field S ought to be in a symmetric product 
of the conjugate of the representation which contains v. If the gauge group is extended 
to contain additional non-abelian factors, such a symmetric product representation will 
generally be large, containing fields which are not neutral under the Standard Model 
gauge group. This jeopardizes asymptotic freedom of the gauge couplings, and in any case 
forces us to view the apparent unification of gauge couplings observed at LEP as merely 
a perverse accident. Furthermore, symmetric product representations are quite difficult 
to obtain in string models. Even worse, the positive gaugino-loop radiative corrections to 
soft scalar masses are proportional to the quadratic Casimir invariant, and thus tend to 
be large for symmetric product representations containing S and S. This effect seems to 
strongly disfavor VEVs for such scalars, although sufficiently large term(s) analogous to 
the A term in (2.2) might overcome it. 

Therefore, we will consider here only the next-to-smallest gauge group containing 
gauged U(1) B -l, namely SU(3) C x SU(2) L x SU(2) R x U(1) B _ L . Under this gauge 
group, the chiral superfields transform as^ 

Q ~ (3, 2, 1,1) L ~ (1, 2, 1, -1) H u , H d ~ (1, 2, 2, 0) 

u,d~ (3, 1,2, -J) e,z/~ (1,1,2,1) 

S -(1,1,3,-2) £-(1,1,3,2). 

The lowest-order non-renormalizable superpotential for the S, S degrees of freedom con- 
tains two independent terms in general: 

W = -^-Tr[SS\ 2 + -^Tr[S 2 ] Tr[S 2 ] . (4.1) 
2M 1 J AM 1 J L J v ; 

Here we use a notation in which SU (2)^> triplets are given by traceless 2x2 matrices, 

explicitly 

S~/y/2 S° \ ( S + /V2 S ++ 

S— S-/V2J ' { S° -S + /V2 / 

with the superscripts indicating the electric charge. The soft breaking terms are given by 
^soft = m 2 Tt[S^S] + m 2 Tr[^£] - {j^^[SS\ 2 + ^Tr[S 2 ] Tr[£ 2 ] + c.c)j . (4.2) 



t This is not a left-right symmetric model, given our (minimal) choice of particle content. 
Similar models are considered in [18], but our treatment will be somewhat different. 
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By a suitable phase rotation, we take the parameter A\ to be real and positive, while the 
phases of Ai, A2, A 2 are arbitrary. There is a possible minimum of the full scalar potential 
for the neutral scalar components of S and S with VEVs in a nearly D-flat direction: 

M 



(S 



0\2 



(a.! + \JA\ - 6|Ai| 2 (m 2 + m 2 )^j , (4.3) 



6|Ai| 2 

with the deviation from D-flatness given by 

(S ) 2 - (S ) 2 * (m 2 - m 2 )/(3g 2 B _ L + 2g\) . 
This minimum is stable against local perturbations provided that 



Af-6|Ai| 2 (m 2 + m 2 ) > 0, 



|Ai + A 2 |V + 



9 9 

m + m 



> 



s 2 \A x + A 2 - 2A|(Ai + X 2 )s\ 2 + 4g R A l 



1/2 



(4.4) 

(4.5) 
(4.6) 



where s = (S°) 2 /M and A 2 = (S°) 2 — (S°) 2 define two convenient parameters of order my/. 
These stability conditions are satisfied in a non- vanishing region of the parameter space. In 
a smaller, but still non- vanishing, region of parameter space this is also a global minimum 
of the potential. (However, it is not clearly relevant to require that the desired minimum 
be global, since the lifetime of the false vacuum might be many orders of magnitude longer 
than the age of the universe.) The VEVs break the gauge symmetry according to 

SU(2) R x U(1) B -l -»■ U(1) Y x (matter parity) . 

A viable spectrum of neutrino masses can arise just as before via the seesaw mechanism, 
by taking the obvious extension of (1.2). The \i term of the MSSM can also arise in 
a way exactly analogous to that discussed in the model of the previous section, from a 
non-renormalizable superpotential term proportional to H u H ( [Yx\SS\. 

The resulting spectrum contains intermediate-scale states consisting of a pair of charge 
±1 massive vector supermultiplets with mass y/2gn(S) and a neutral massive vector su- 
permultiplet of mass yj §g 2 B _ L + (S) ■ The remaining uneaten components of S and S 
obtain electroweak-scale masses, and consist of the pair of charge ±2 chiral supermultiplets 
S , S ++ and one neutral chiral supermultiplet $. The components of $ get masses just 
given by eqs. (2.10)-(2.11) with A — > Ai and A — > A\; and all of the same comments apply 



as before. The fermionic components of S and S 
their light scalar partners have squared masses 



T7+ + 



obtain masses |Ai + A2 1 s, while 



|Ai+A 2 |V + 



9 9 

m + m 



± 



s 2 \A± +A 2 - 2Af(Ai + X 2 )s\ 2 + 4g A R A l 
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1/2 



(4.7) 



These masses suffer renormalization between the scale of spontaneous symmetry breaking 
and the electroweak scale, since S , S are charged under the MSSM gauge group. 

The most striking prediction of this model is therefore the presence of an exotic vector- 
like pair of chiral supermultiplets of electric charge ±2 which may well be accessible to 
future collider experiments. These particles will have unsuppressed two-body decays into 
pairs of like-sign leptons, because of the superpotential interaction W D ys ee -S which 
derives from the analog of (1.2). This should yield a striking experimental signature. 
This feature is shared by left-right symmetric models with symmetry breaking near the 
electroweak scale [33]. In the present case, the lightness of these exotic states is due to 
the lack of renormalizable mass couplings in the underlying superpotential. One can check 
that the presence of these exotic states does not cause the gauge couplings to blow up 
below the Planck scale; however, they do completely modify the running. In this model, 
gauge coupling unification in the usual sense would require additional fields not considered 
here, and unlike in the model of section 2, the LEP observation of apparent unification 
would have to be viewed as entirely accidental. 



5. Conclusion 



In this paper we have analyzed what might be called the minimal supersymmetric stan- 
dard model with automatic .R-parity conservation. We do not include renormalizable tree- 
level mass terms in the superpotential. Instead, gauge symmetry breaking arises because 
of the interplay between non-renormalizable interactions and soft supersymmetry-breaking 
interactions. We then found that it is important to take into account dimensionless but 
"soft" supersymmetry-breaking couplings in this analysis, which can play a crucial role in 
the spontaneous symmetry breaking; indeed, it may not be able to understand the sym- 
metry breaking mechanism without them. These models have several attractive features. 
First, the \i term of the MSSM is naturally generated by a mechanism familiar from [27]. 
Second, the masses of neutrinos are determined by an intermediate-scale seesaw mechanism 
and so may be phenomenologically interesting. The minimal version of the model in sec- 
tion 2 also has the nice property that the apparent unification of gauge couplings observed 
at LEP can still be considered non-accidental. In this model we found that, assuming 
supergravity-inspired boundary conditions on the soft terms, a discernible imprint may 
be left on the spectrum of MSSM sparticles. In particular, the masses of the left-handed 
sleptons are further increased over those of the right-handed sleptons. 
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One of the interesting consequences of this class of models is the existence of a su- 
permultiplet (<&) of neutral particles with electroweak scale masses and only very weak 
couplings to MSSM particles. The largest couplings of these particles to MSSM states 
are suppressed by at least (m^/M) 1 / 2 , and arise from the same non-renormalizable in- 
teraction which induces the ji term of the MSSM. These particles should therefore have 
relatively long (perhaps microsecond or nanosecond) lifetimes. While this may provide for 
interesting cosmological consequences, the weak couplings of these particles means that 
they cannot play a role in collider experiments. In the next-to-minimal model, we found 
that the symmetry breaking mechanism also predicts a pair of exotic chiral supermultiplets 
of particles with electric charge ±2 and electroweak-scale masses. 

It is possible that -R-parity conservation cannot be "explained" , but should simply be 
taken as a law of nature. It is also possible that an explanation exists, but lies only on the 
far side of the Planck or string scale. However, it is gratifying that one can construct field 
theory models which are consistent with all known observations, and in which .R-parity 
conservation has its origin in terms of a deeper gauge principle. This may be taken as one 
of many clues to the nature of physics at very high energies. 

I am grateful to James D. Wells for many helpful discussions. This work was supported 
in part by the U.S. Department of Energy. 
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